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The inclusive cross section for production of isolated photons has been measured in pp collisions 
at y/s — 630 GeV with the D0 detector at the Fermilab Tevatron Collider. The photons span a 
transverse energy (Et) range from 7-49 GeV and have pseudorapidity |r?| < 2.5. This measurement 
is combined with the previous D0 result at y/s — 1800 GeV to form a ratio of the cross sections. 
Comparison of next-to-leading-order QCD with the measured cross section at 630 GeV and the ratio 
of cross sections show satisfactory agreement in most of the Et range. 
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Within the framework of Quantum Chromodynamics 
(QCD), isolated single photons are direct photons: pro- 
duced from the primary parton-parton interactions. Be- 
cause the dominant production mechanism for photons 
of modest transverse energy (Et) at the Fermilab Teva- 
tron is gluon Compton scattering (qg — > 7(7), the cross 
section for direct-photon production is sensitive to the 
gluon distribution in the proton jjj. A measurement of 
the final state photons provides a probe of QCD with- 
out additional complications from fragmentation and jet 
identification, providing a powerful and effective means 
for studying the constituents of hadronic matter. 

Previous experiments, at center-of-mass energies of 
both 630 GeV § and 1800 GeV §@, have reported pho- 
ton production in excess of next-to- leading-order (NLO) 
QCD predictions at low transverse energies (Ej, < 30 
GeV). This disagreement with data could result from 
gluon radiation not included in NLO calculations || or 
because the parton distributions are not well known 

In this Letter, we present a measurement of the iso- 
lated photon cross section in pp collisions for photons in 
two pseudorapidity regions, \r]\ < 0.9 and 1.6 < |ry| < 2.5, 
where r\ = — In tan | and 9 is the polar angle with respect 
to the proton beam. We compare the production cross 
section at \fs — 630 GeV with the previously published 
D0 results at -y/s = 1800 GeV f§. A ratio of the cross 
sections at different energies reduces systematic uncer- 
tainties and minimizes the sensitivity to the choice of 
parton distribution functions (PDF) because the mea- 
surements at both energies use the same detector and 
the same analysis method. 

The cross section measurement at 630 GeV uses a sam- 
ple of 520 nb _1 of data recorded in 1995 [jjj with the D0 
detector at the Fermilab Tevatron Bj . The analysis uses 
the uranium/liquid argon calorimeter to identify electro- 
magnetic (EM) showers, and the drift chambers in front 
of the calorimeter to differentiate photon showers from 
electron showers. The EM calorimeter provides full az- 
imuthal ((f>) coverage, and consists of a central cryostat 
(CC) with |7y| < 1.1, and two forward cryostats (EC) with 
1-4 ^ \r)\ ^ 4.0. The EM calorimeter is divided into four 
longitudinal layers, EM1-EM4, of approximately 2, 2, 7, 
and 10 radiation lengths, respectively. The EM energy 
resolution in the central and forward calorimeter is given 
by a E /E = {15%/v/£(GeV)} © 0.3%. 

Photons interacting in the calorimeter are detected us- 
ing a three-level triggering system. The first level con- 
sists of scintillation counters near the beam pipe, which 
detect inelastic pp collisions. The second level requires 
a minimum energy deposition in a A(f> x A77 = 0.2 x 0.2 
trigger tower, with thresholds of 2.0, 3.0, and 7.0 GeV. 
In the final step, calorimeter clusters are formed with 
corresponding thresholds of 4.5, 8.0, and 14.0 GeV. The 
trigger efficiency is determined for the 14.0 and 8.0 GeV 
thresholds by taking the ratio of events passing each trig- 
ger criteria to those passing the 8.0 and 4.5 GeV criteria, 



respectively, in an energy regime where the lower thresh- 
old trigger is 100% efficient. Monte Carlo studies of the 
trigger algorithms show agreement with the data for the 
two higher energy triggers, and are used to determine 
the trigger efficiency for the 4.5 GeV trigger. Trigger ef- 
ficiencies are typically about 20% at the nominal energy 
threshold and rise to almost 100% a few GeV above the 
threshold value. Consequently, photon candidates are 
accepted only for transverse energies of at least 7.35, 10, 
and 16 GeV for the three triggers, respectively. 

Photon candidates are identified as energy clusters lo- 
cated well within the pseudorapidity boundaries of the 
central calorimeter or the forward calorimeter, and, in the 
central calorimeter, located at least 1.6 cm from the az- 
imuthal section boundaries. The event vertex position is 
required to be within 50 cm of the center of the detector. 
The resulting geometric acceptance is A = 0.622 ± 0.007 
(0.787 ± 0.007) in the central (forward) region. Candi- 
dates must pass a series of selection criteria || , that iden- 
tify the energy cluster as an electromagnetic shower. The 
total transverse energy near any candidate cluster must 
satisfy an isolation requirement EJ^~ 0A — E^- 02 < 2.0 
GeV, where 1Z = J (A77) 2 + (A<f>) 2 is the distance from 
the cluster center. The combined selection and isolation 
efficiency, e s , is estimated as a function of E^, from a 
GEANT-based Monte Carlo simulation of the D0 detec- 
tor. We find e s ~ 60% (75%) in the CC (EC) at 8.0 GeV 
and e s - 88% (90%) above 20 GeV. To minimize back- 
ground from electrons, photon candidates are rejected if 
any tracks in the drift chamber extrapolate to within a 
road of width A</> x A9 = 0.2 x 0.2 defined by the angle 
subtended by the candidate photon cluster and the initial 
interaction vertex. The total charged tracking efficiency 
is estimated from Z — > e + e~ events to be 0.858 ± 0.013 
(0.593 ± 0.079) in the central (forward) region. 

The predominant background to direct photon pro- 
duction arises from the decay of ir° or 77 mesons to two 
photons. The fraction of direct photons is determined 
from the energy (Ei) deposited in the innermost longi- 
tudinal section of the calorimeter, EMI. Photons have 
a small probability of showering in the material in front 
of the calorimeter and, thus, tend to deposit little en- 
ergy in EMI. Sensitivity to the amount of EMI en- 
ergy can be used to distinguish multiple photon back- 
ground from a single photon signal. We use the function 
f{Ei) = log 10 [l + log 10 {l -I- i^i(GeV)}] as our discrimi- 
nant to determine the single photon purity. The expected 
distributions of this function for signal and background 
are found from events simulated with the PYTHIA Monte 
Carlo |j] and overlaid with data acquired using a random 
trigger to model noise, pileup, and multiple pp interac- 
tions. Three categories of fully simulated events are gen- 
erated: those containing photons, and background events 
with and without charged tracks pointing from the in- 
teraction vertex to the EM cluster. The two different 
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background samples are generated so that charged and 
neutral background fractions can be separately fit to the 
data, thus minimizing uncertainties from the tracking ef- 
ficiency and from the model used for jet fragmentation. A 
systematic uncertainty in modeling jet fragmentation is 
estimated by varying the multiplicity of neutral mesons in 
the core of PYTHIA jets by ±10%. The detector response 
is modeled using a detailed GEANT simulation with the 
energy response in EMI calibrated to match the data 
from W — ► ev events. 

The same criteria used to select photon candidates in 
the data are applied to the Monte Carlo events. The dis- 
tribution of / from the data is fitted to a normalized lin- 
ear combination of Monte Carlo photons and background 
with and without charged tracks in the road pointing 
back to the interaction vertex. The fit is performed in 
different EZ regions using the CERNLIB fitting package 
hmcmll |1(| , with the fractions of signal and background 
constrained to be between 0.0 and 1.0. The purity is de- 
fined as the fraction of Monte Carlo photons in the nor- 
malized fitted distribution. A representative fit is shown 
in Fig. U and the photon purity as a function of E^ is 
plotted in Fig. §. 
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FIG. 1. Distribution of the discriminant, /, for determining 
photon purity, where E\ is in units of GeV. Points with error 
bars indicate data. Broken lines indicate simulated distribu- 
tions of (a) single photons, and jet background (b) without 
and (c) with charged tracks. The solid line depicts a fit sum 
of all three distributions. 

The final cross sections d 2 'a / 'dE^drj , after applying ef- 
ficiency and purity corrections, are shown in Fig. |^ and 
tabulated in Table |. The error bars show all uncorre- 
cted uncertainties, which include the statistical uncer- 
tainty, and uncertainties from selection criteria, trigger 
efficiency, and the fitted photon purity. The contribu- 
tion from the fit to photon purity is the largest source 
of uncorrelated uncertainty. The correlated uncertainty 
consists of the uncertainties in luminosity, tracking ef- 
ficiency, geometric acceptance, calorimeter energy scale, 
and the largest contribution, that from the fragmentation 
model. 
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FIG. 2. The photon purity as a function of E^, for central 
and forward photons. The error bars indicate the uncertainty 
in the fit purity and are larger than errors derived from sta- 
tistical analysis alone. 
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FIG. 3. The cross section for production of isolated pho- 
tons in central and forward regions at ^J~s — 630 GeV. The 
error bars show the total uncorrelated error, and the curves 
show cross sections predicted from NLO QCD. 
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TABLE I. The measured and predicted isolated photon 
production cross section at yHi = 630 GeV. The value for 
the column labeled "Plotted _BJ" is determined according to 
Ref. |l^] . The columns labeled Sau and Sac are, respectively, 
the uncorrelated and correlated uncertainties. 
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The results are compared with NLO QCD calculations 
using CTEQ5M parton distributions jy]], with renormal- 
ization and factorization scales fj,R = /i^ = _E™ ax , where 
£"™ ax is the maximum photon transverse energy in the 
event. Figure ^ compares the data and theory. A co- 
variance matrix x 2 provides a measure of the probability 
that the theory describes the data. A complete covari- 
ance matrix, composed of correlated and uncorrelated 
uncertainties, is determined and the theoretical cross sec- 
tion is compared to the data with a % 2 value of 11 (4.6) 
for 7 degrees of freedom in the CC (EC) region. This 
gives a standard % 2 probability that the theory is con- 
sistent with the data at 12% (71%) probability in the 
CC (EC) regions. Deviations between theory and data 
are largest at low Ej, in the central region. These re- 
sults are in qualitative agreement with those previously 
published at y/s = 1800 GeV, where the theory is lower 
than the data at low E T (wlO-40 GeV) in the CC, but 
is consistent with the data over all E T in the EC ||. 
Using different PDFs changed the cross section by less 
than 5% |u|. Setting scales to fx R = fj, F = 2.0£™ ax or 
/itR = [If = 0.5i?™ ax changed the cross section by about 
20%, as shown in Fig. [|. 
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FIG. 4. Comparison of the measured cross section for pro- 
duction of isolated photons at yfs = 630 GeV with NLO QCD 
using CTEQ5M parton distribution functions. The error bars 
indicate the uncorrelated uncertainty and the shaded bands 
indicate the correlated uncertainty 

In the simple parton model, the dimensionless cross 



section E T 



E'- 



as a function of xt — 2 -j=-, is in- 



dependent of \fs. Although deviations from such naive 
scaling are expected, the dimensionless framework pro- 
vides a useful context for comparison with QCD. The 
experimental dimensionless cross section, averaged over 
azimuth, becomes op = ^ • d 2 a / 'dExdr/. The ratio 
a D {y/s = 630 GeV)/a D (y/s = 1800 GeV) is determined 
by combining the cross section reported in this Letter 
with the D0 measurement at s/s = 1800 GeV 



The ratio is shown as a function of xt in Fig. and 
Table |l| together with the NLO QCD prediction. 
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FIG. 5. The ratio of the dimensionless cross sections, 
(T D (s/s = 630 GeV)/a D (y r s = 1800 GeV). The error bars 
indicate the uncorrelated uncertainty and the shaded bands 
indicate the correlated uncertainty. 

Comparison of the theoretical cross section ratio to 
the data, using the complete covariance matrix, gives a 
X 2 value of 6.5 (3.0) for 7 degrees of freedom in the CC 
(EC), which corresponds to a standard x 2 probability of 
49% (89%) in the CC (EC) region. Although the low- 
est xt points are systematically higher than NLO QCD 
predictions in both the CC and EC regions, the devia- 
tions are not significant in light of our combined statis- 
tical and systematic uncertainties, and there exists good 
agreement between the measured ratio and theory. 

We have measured the production cross section for iso- 
lated photons in pp collisions at s/s = 630 GeV and com- 
pared this cross section with that measured at s/s = 1800 
GeV. The measurement is higher than the theoretical 
prediction at low Et in the central rapidity region but 
agrees at all other Et and in the forward rapidity region. 
The difference between data and theory is less significant 
for the ratio of cross sections, and the theory is consistent 
with the data over all Et- 
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that their calculated cross sections are consistent over the 
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TABLE II. The measured ratio and NLO QCD prediction 
for the dimensionless cross section at y^s = 630 GeV to that 
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